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ABSTRACT 

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. This is the first of two papers where we 
present the results of a thorough temporal and spectral analysis of all the data available to us, focusing in particular 
on the flare of April 1998, which was simultaneously observed also at TeV energies. Here we focus on the time 
analysis, while the spectral analysis and physical interpretation are presented in the companion paper. The detailed 
study of the flare in different energy bands reveals very important new results: i) hard photons lag the soft ones 
by 2-3 ks - a behavior opposite to what is normally found in high energy peak BL Lacs X-ray spectra; ii) the 
flare light curve is symmetric in the softest X-ray band, while it becomes increasingly asymmetric at higher 
energies, with the decay being progressively slower than the rise; iii) the flux decay of the flare can be intrinsically 
achromatic if a stationary underlying emission component is present. The temporal and spectral information 
obtained challenge the simplest models currently adopted for the (synchrotron) emission and most importantly 
provide clues on the particle acceleration process. 

Subject headings: galaxies: active — BL Lacertae objects: general — BL Lacertae objects: individual (Mkn 
421) — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

Blazars are radio-loud AGNs characterized by strong vari- 
ability, large and variable polarization, and high luminosity. Ra- 
dio spectra smoothly join the infrared-optical-UV ones. These 
properties are successfully interpreted in terms of synchrotron 
radiation produced in relativistic jets and beamed into our di- 
rection due to plasma moving relativistic ally close to the line 
of sight (e.g. Urry & Padovani 1995). Many blazars are also 
strong and variable sources of GeV 7-rays, and in a few ob- 
jects the spectrum extends up to TeV energies. The hard X- to 
7-ray radiation forms a separate spectral component, with the 
luminosity peak located in the MeV-TeV range. 

The emission up to X-rays is thought to be due to syn- 
chrotron radiation from high energy electrons in the jet, while it 
is likely that 7-rays derive from the same electrons via inverse 
Compton (IC) scattering of soft (IR-UV) photons -synchrotron 
or ambient soft photons (e.g. Sikora, Begelman & Rees 1994, 
Ghisellini et al. 1998). 

The contributions of these two mechanisms characterize the 
average blazar spectral energy distribution (SED), which typi- 
cally shows two broad peaks in a vF v representation (e.g. von 
Montigny et al. 1995; Sambruna, Maraschi & Urry 1996; Fos- 
sati et al. 1998a): the energies at which the peaks occur and 
their relative intensity provide a powerful diagnostic tool to in- 
vestigate the properties of the emitting plasma, such as electron 
energies and magnetic field (e.g. Ghisellini et al. 1998). More- 
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over variability studies, both of single band and of simultane- 
ous multifrequencies data, constitute the most effective means 
to constrain the emission mechanisms at work in these sources 
as well as the geometry and modality of the energy dissipation. 
The quality and amount of X-ray data on the brightest sources 
start to allow us to perform a thorough temporal analysis as 
function of energy and determine the spectral evolution with 
good temporal resolution. 

In X-ray bright BL Lacs (HBL, from High-energy-peak-BL 
Lacs, Padovani & Giommi 1995) the synchrotron maximum 
(usually) occurs in the soft-X-ray band, and the inverse Comp- 
ton emission extends in some cases to the TeV band where - 
thanks to ground based Cherenkov telescopes - four sources 
have been detected up to now: Mkn 421 (Punch et al. 1992), 
Mkn 501 (Quinn et al. 1996), 1ES 2344+514 (Catanese et al. 
1998), and PKS 2155-304 (Chadwick et al. 1998). If the inter- 
pretation of the SED properties in terms of synchrotron and IC 
radiation is correct, a correlation between the X-ray and TeV 
emission is expected. 

Mkn 421 (z = 0.031) is the brightest BL Lac object at X- 
ray and UV wavelengths and the first extragalactic source dis- 
covered at TeV energies, where dramatic variability has been 
observed with doubling times as short as 15 minutes (Gaidos 
et al. 1996). As such it was repeatedly observed with X-ray 
satellites, including BeppoSAX. Remarkable X-ray variability 
correlated with strong activity at TeV energies has been found 
on different occasions (Macomb et al. 1995, 1996, Takahashi et 
al. 1996, Fossati et al. 1998b, Maraschi et al. 1999). In partic- 
ular, the 1998 BeppoSAX data presented here were simultane- 
ous with a large TeV flare detected by the Whipple Observatory 
(Maraschi et al. 1999). 

This paper is the first of two, which present the results of a 
uniform, detailed spectral and temporal analysis of BeppoSAX 
observations of Mkn 421 performed during 1997 and 1998. 
Here we focus on the data reduction and the timing analysis, 
and also discuss the results on the spectral variability derived 
from the different properties of the flux variations in different 
energy bands. 
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Evolution of Mkn 421: I. Timing Analysis 



The paper is organized as follows. We briefly summarize the 
characteristics of BeppoSAX (§2), and introduce the observa- 
tions studied (§3). We then address the temporal analysis of 
the variability, considering several energy bands and compar- 
ing the light curve features by means of a few simple estima- 
tors for the 1997 and 1998 observations (§4). The remarkable 
flare observed in 1998 is the object of a further deeper analysis, 
reported in Section §5, focused on timescales and time lags. 
Section 6 contains a summary of the results of the temporal 
analysis, preparing the ground for the comprehensive discus- 
sion presented in Paper II (Fossati et al. 2000). There they are 
considered together with the results of the spectral analysis and 
thus used to constrain a scenario able to interpret the complex 
spectral and temporal findings. 

2. BEPPOSAX OVERVIEW 

For an exhaustive description of the Italian/Dutch BeppoSAX 
mission we refer to Boella et al. (1997) and references 
therein. The narrow field coaligned instrumentation (NFI) on 
BeppoSAX consists of a Low Energy Concentrator Spectrom- 
eter (LECS), three Medium Energy Concentrator Spectrome- 
ters (MECS), a High Pressure Gas Scintillation Proportional 
Counter (HPGSPC), and a Phoswich Detector System (PDS). 
The LECS and MECS have imaging capabilities in the 0.1- 
10 keV and 1.3-10 keV energy band, respectively, with en- 
ergy resolution of 8% at 6 keV. At the same energy, the an- 
gular resolution is about 1.2 arcmin (Half Power Radius). In 
the overlapping energy range the MECS effective area (150 
cm 2 ) is ~ 3 times that of the LECS. Furthermore the expo- 
sure time for the LECS is limited by stronger operational con- 
straints to avoid UV light contamination through the entrance 
window (LECS instrument is operated during Earth dark time 
only). The HPGSPC covers the range 4-120 keV, and the PDS 
the range 13-300 keV. In the overlapping interval the HPGSPC 
has a better energy resolution than the PDS, but it is less sensi- 
tive of both PDS and MECS. Therefore, HPGSPC data will not 
be discussed in this paper. 

The present analysis is based on the SAXDAS linearized 
event files for the LECS and the three MECS experiments, to- 
gether with appropriate background event files, as produced at 
the BeppoSAX Science Data Center (REV 0.2, 1.1 AND 2.0). 
The PDS data reduction was performed using the XAS software 
(Chiappetti & Dal Fiume 1997) according to the procedure de- 
scribed in Chiappetti et al. (1999). 

3. OBSERVATIONS 

Mkn 421 has been observed by BeppoSAX in the springs of 
1997 and 1998. The journal of observations is given in Table 1. 

3.1. 1997 

The 1997 observation comprised several pointings spanning 
the interval between April 29 th and May 7 th . MECS data are not 
available for May 7 th because of the failure of the detector unit 
1 on May 6 th . In this paper we will not consider the LECS data 
of this last day of the 1997 campaign, because unfortunately 
LECS data alone do not provide useful spectral and variabil- 
ity information. The net exposure time (excluding the May 7 th 
data) was ~ 52 and ~ 1 17 ks for LECS and MECS respectively, 
while the on-source time coverage computed as the sum of the 
observations between each T start and T stop is of about 58 hr. 

Results on the first half of the 1997 campaign (April 29 th 
to May 1 st ) have been presented by Guainazzi et al. (1999), 



where the details and motivation of the observations are given. 
We re-analyzed those data along with the new ones, applying 
the same techniques, in order to obtain a homogeneous set of 
results, necessary for a direct comparison. 

3.2. 1998 

In 1998 BeppoSAX observed Mkn 421 as part of a long 
monitoring campaign involving BeppoSAX, ASCA (Takahashi, 
Madejski, & Kubo 1999), RossiXTE (Madejski et al., in prepa- 
ration), and coverage from ground based TeV observatories 
(Maraschi et al. 1999). The BeppoSAX observation comprised 
two distinct long pointings, started respectively on April 21 st 
and 23 rd . The total net exposure time for MECS telescopes has 
been of 64 ks and the LECS one adds up to 50 ks. Hereinafter 
we will refer to the two 1998 pointings simply as April 21 st and 
April 23 rd , leaving out the year. 

The actual on-source time was 55.5 hr. Unfortunately celes- 
tial and satellite constraints in late April 1998 were such that 
in each orbit there are two intervals during which the recon- 
structed attitude is undefined. One occurs when the source is 
occulted by the dark Earth, but the second longer one occurs 
during source visibility periods. Thus there are about 19 min 
per orbit when the NFIs are pointing at the source, but there 
is a gap in attitude reconstruction despite the high Earth eleva- 
tion angle. During these intervals the source position is drift- 
ing along the ±X detector axis, from a position in the center 
to about 12 arcmin off, passing the strongback window sup- 
port. These intervals are excluded by the LECS/MECS event 
file generation. Using the XAS software it is possible to accu- 
mulate MECS images also during such intervals, showing that 
the source is actually drifting. The drift causes a strong energy 
dependent modulation, clearly visible in light curves accumu- 
lated with XAS, and which would be very difficult to recover, 
as requires to build a response matrix integrated from differ- 
ent position dependent matrices weighted according to the time 
spent in each position. We then did not try to recover these data 
for the analysis. 

In the case of PDS the source remains inside the collimator 
flat-top during the drift, and therefore we used for the accu- 
mulation all events taken when the Earth elevation angle was 
greater than 3 degrees (inclusive of attitude gaps). 

4. TEMPORAL ANALYSIS 

Light curves have been accumulated for different energy 
bands. The input photon lists were extracted from the full 
dataset selecting events in a circular region centered on the posi- 
tion of the point source. The extraction radii are 8 and 6 arcmin 
for LECS and MECS respectively, chosen to be large enough 
to collect most of the photons in the whole energy range 9 : the 
used values ensure that < 5 % of the photons are missed, at all 
energies. 

For MECS data we considered the merged photon list of all 
the available MECS units, which were 3 for the 1997 observa- 
tions and only 2 (namely MECS units 2 and 3) in 1998. 

The expected background contribution is less than a fraction 
of a percent (typically of the order of 0.2-0.5 %) and therefore 
the light curves have not been background subtracted. 

In order to maximize the information on the spectral/temporal 
behavior, the choice of energy bands shall be such that they are 

9 For bright and soft sources - this matters especially because the LECS 
point spread function (PSF) gets rapidly broad below the Carbon edge, i.e. E < 
0.3 keV - it is recommended to select quite large values. 
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Table 1 
Journal of Observations 



Date (UTC) Net Exposure Time (ks) BeppoSAX archive # 



Start 


Stop 


LECS 


MECS 


PDS a 




1997/04/29:04:02 


1997/04/29:14:42 


11.6 


21.8 


18.4 


50032001 


1997/04/30:03:19 


1997/04/30:14:42 


11.4 


24.1 


23.0 


50032002 


1997/05/01:03:17 


1997/05/01:14:42 


11.2 


23.8 


22.8 


50032003 


1997/05/02:04:10 


1997/05/02:09:41 


4.4 


11.4 


11.0 


50016002 


1997/05/03:03:24 


1997/05/03:09:41 


4.3 


11.7 


11.2 


50016003 


1997/05/04:03:25 


1997/05/04:09:45 


4.9 


12.2 


11.6 


50016004 


1997/05/05:03:32 


1997/05/05:09:45 


5.0 


11.9 


11.4 


50016005 


1997/05/07:04:47 


1997/05/07:10:27 


6.0 




13.1 


50016006 b 


1998/04/21:01:52 


1998/04/22:03:13 


23.6 


29.6 


42.2 


50686002 


1998/04/23:00:27 


1998/04/24:06:37 


27.2 


34.7 


50.4 


50686001 



a Net exposure time over half of the total area, because of the rocking mode in which the instrument is 
operated. 

b On May 7 th 1997 MECS 1 had a fatal failure. Only LECS and PDS data are available for that day, and 
they are not discussed in the paper. 



as independent as possible. This should take into account the 
instrumental efficiencies and, at a certain level, the details of the 
spectral shape. The main trade-off is the number of photons in 
each band, that should be large enough to obtain statistically 
meaningful results. Taking into account the main features of 
the LECS and MECS effective areas (e.g. the Carbon edge at 
0.29 keV which provides with independent energy bands be- 
low and above this energy) and the steepness of the Mkn 421 
spectrum (above a few keV the available photons quickly be- 
come insufficient), we determined the following bands: 0.1- 
0.5 keV and 0.5-2 keV for LECS data, 2-3 and 4-6 keV for 
MECS data. Their respective "barycentric" energies 10 are 0.26, 
1.16, 2.35, 4.76 keV, respectively, thus providing a factor ~ 20 
leverage, useful to test the energy dependence of the variability 
characteristics. 

4.1. 1997 

The four resulting light curves for 1997 are presented in Fig- 
ure 1 . The source shows a high degree of flux variability, with 
possibly a major flare between the third and fourth pointings. 
The vertical scale is logarithmic allowing a direct comparison 
of the amplitude of variations at different energies. As antic- 
ipated we are not going to consider further the LECS data of 
May 7 th , however for sake of completeness we report that the 
count rate was of 0.37±0.01 and 1.06±0.04 cts/s for the 0.1- 
0.5 and 0.5-2.0 keV energy bands respectively, comparable to 
the level measured on May 5 th . 

The comparison of data with the overlaid reference grid 
shows that the amplitude of variability is larger at higher en- 
ergies, as commonly observed in blazars at energies above the 
SED peak (e.g. Ulrich, Maraschi & Urry 1997). We will dis- 
cuss this issue more quantitatively in section §4.3. 

4.2. 1998 

The 1998 light curves are shown in Figure 2, (see the caption 
for details about the axis scales). The overall appearance is 

10 These have been computed directly from the count spectra and thus they 
are not model dependent. Moreover, despite of the observed spectral variability, 
their values do not change more than a few eV. 



Date of 1997 




2x10 s 4xl0 5 

Time [s from T = 10567.0 (TJD)] 

FIG. 1 . — Light curves in four energy bands for the 1997 observations. They 
are shown in order of decreasing energy from top to bottom: 4-6 keV (mul- 
tiplied by a factor 100), 2-3 keV (multiplied by a factor 10), 0.5-2 keV, 0.1- 
0.5 keV. The binning time is 1500 s. The reference time TJD=10567.0 corre- 
sponds to 1997-April-29:()0:()0 UT. The horizontal bars shown at the bottom 
of the plot mark the time intervals considered for the analysis of the variability 
properties in §4.3: subset-1 to -3 from left to right. 



somewhat different, being dominated by a single isolated flare 
at the beginning of the campaign. The variability amplitude is 
similar to the 1997 one in each corresponding energy band (see 
also §4.3). 

The most striking and important result of the campaign is 
that in correspondence with the X-ray flare of April 21 st a 
sharp TeV flare was detected by the Whipple Cherenkov Tele- 
scope, with amplitude of a factor 4 and a halving time of about 
2 hr (Maraschi et al. 1999). In Figure 3 the Whipple TeV 
(E>2 TeV) light curve is shown together with the BeppoSAX 
ones: LECS (0.1-0.5 keV), MECS (4-6 keV) and PDS (12- 
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Date of 1998 

] April 21 ] April 22 ] April 23 ] April 24 J 

I I I I I I 




2xl0 5 4xl0 5 

Time [s from T = 10924.0 (TJD)] 

FIG. 2. — Light curves in four energy bands for the 1998 observations. They 
are shown in order of decreasing energy from top to bottom: 4-6 keV (mul- 
tiplied by a factor 100), 2-3 keV (multiplied by a factor 10), 0.5-2 keV, 0.1- 
0.5 keV. The binning time is 1500 s. The reference time TJD=10924.0 cor- 
responds to 1998-April-21:00:00 UT. In order to allow an easier comparison 
with the 1997 light curves, the X scale has been expanded to the same length 
of that of Figure 1, and also the Y-axis dynamic range is the same. 



Date of 1998 

April 21 April 22 April 23 
500 I — i 1 1 




10 s 2xl0 6 3xl0 5 

Time [s from T = 10924.0 (TJD)] 



FIG. 3. — Light curves at TeV and X-ray energies, during the 1998 campaign. 
They are shown in order of decreasing energy from top to bottom: Whipple E> 
2 TeV, (multiplied by a factor 500), PDS 12-26 keV (one point per orbit, multi- 
plied by a factor 50), MECS 4-6 keV (multiplied by a factor 5) and LECS 0. 1- 
0.5 keV (both with 1500 s bins). The count rate units are cts/s for BeppoSAX 
data, and cts/min for Whipple data. The thicker dashed grid marks the days, 
while the thinner ones have a 20 ks spacing. TJD= 10924.0 corresponds to 
1998-April-21:00:00 UT. 



26 keV) instruments. The peaks in the 0.1-0.5 keV, 4-6 keV 
and 2 TeV light curves are simultaneous within one hour. Note 
however that the TeV variation appears to be both larger and 
faster than the X-ray one. The 12-26 keV light curve shows a 
broader and slightly delayed peak with respect to lower energy 
X-rays. However the very limited statistics does not allow us 
to better quantify this event. A more detailed account on this 
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FIG. 4. — RossiXTE All Sky Monitor light curve (top) compared with the 
MECS in the same (2-10 keV) band, scaled by a factor 1/20. 



particular result is presented and discussed by Maraschi et al. 
(1999). 

It is also worth noticing that the positive (flux) offset of the 
beginning of April 23 rd with respect to the end of the April 
21 st suggests the presence of a second flare occurring between 
the two pointings. This seems to be also confirmed by the 
RossiXTE All Sky Monitor (ASM) light curve, shown in Fig. 4 
together with the MECS one in the (2-10 keV) nominal work- 
ing energy range of the ASM. Single ASM ^90 s dwells have 
been rebinned in 14400 s (i.e. 4 hr) intervals, weighting the 
contribution of each single dwell on its effective exposure time 
and on the quoted error. The MECS light curve bin is 1000 s. 
Indeed the ASM light curve appears to unveil the presence of 
a second flare occurring in between the two BeppoSAX point- 
ings, with a brightness level similar to that of the detected one. 

On the other hand the few Whipple data points at the time of 
this putative second X-ray flare (at T = (100-120) x 10 3 s, see 
Figure 3) do not indicate any (major) TeV activity: the count 
rate measured by the Whipple telescope was in fact signifi- 
cantly lower than that measured simultaneously with the first 
X-ray flare. 

4.3. Comparison of 1997 and 1998 variability characteristics 

We can characterize the variability in different X-ray bands 
by two commonly used estimators: the fractional root mean 
square (r.m.s.) variability parameter F var , and the minimum 
"doubling/halving time " T s h rt (for definitions and details see 
Appendix A, and Zhang et al. 1999). 

4.3.1. Fractional r.m.s. variability 

For each (0. 1-0.5, 0.5-2, 2-3 and 4-6 keV) band we computed 
F var for the whole 1997 and 1998 datasets which have the same 
on-source coverage (i.e. ~ 55 hr). 

Moreover, for the 1997 we considered three (partially over- 
lapping) subsets of data each spanning a time interval compara- 
ble to the length of the whole 1998 observation, i.e. ~ 300 ks. 
The intervals cover the first 300 ks, a middle section, and the 
last 300 ks, as marked on Fig. 1. 
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FIG. 5. — Plot of the root mean square variability parameter, F va r, versus 
energy. A small energy shift has been applied for sake of clarity. 



For each dataset F var is computed for 6 different light curve 
binning intervals, namely 200, 500, 1000, 1500, 2000 and 2500 
s, and then the weighted average 11 of these values is computed. 
The results are listed in Table 2 and shown in Figure 5. 

If we treat all the datasets as independent and momentarily 
do not consider the subset-1 of the 1997 light curve (empty 
triangles in Fig. 5) on which we will comment later, it appears 
that: 

1 . In terms of root mean square variability the lower energy 
flux is less variable than the higher energy one, as already 
pointed out in section §4.1. This holds independent of 
both the length of the time interval and the state of the 
source. 

2. There is no relation between the brightness level and the 
variability amplitude, since although in 1998 the source 
was at least a factor 2 brighter than in 1997 the corre- 
sponding F var are not larger. 

What is different for the subset-1 of the 1997 light curve ? 

First two caveats on F var . This quantity is basically only sen- 
sitive to the average excursion around the mean flux, and does 
not carry any information about either duty cycle or the actual 
flux excursion. Therefore in order for F var to meaningfully rep- 
resent the variability, the "extrema" of the source variability 
have to be sampled. This is clearly not the case for the first 1997 
sub-light curve (see Figure 1). Furthermore, F var is estimated 
from the comparison of the variances of the light curve and the 
measurements, but while the latter ones probably obey a Gaus- 
sian distribution and are much smaller than the former ones, 
the "probability" for the source to show a significant deviation 
from the average (i.e. the count rate histogram distribution) 
is in most cases highly non-Gaussian. For a typical "flaring" 
light curve this has a "core" at low count rates, resulting from 
the long(er) time spent by the source in between flares, and a 
very extended tail to high(er) rates. The implicit assumption of 
Gaussian variability can thus affect the estimate of F var (as for 
large amplitude and low duty cycle variability). 

"To the F var values obtained for each different time binning, we assigned as 
weight = 1 /ct 2 , where a is their respective one sigma uncertainty. 



F var is thus a particularly poor indicator especially for obser- 
vations with a window function like the one of 1997, for which 
the source temporal covering fraction is of the same order of 
the single outburst, making likely to miss the flare peak, and 
the total span of the observation does not significantly exceed 
the frequency of peaks. We therefore believe this affects the 
results relative to the first subsample of 1997. 

In order to obtain a variability estimator more sensitive to 
the dynamic range towards higher count rates (flaring) we also 
considered a modified definition of it, F var mec i: to study only 
the width (amplitude) of the brighter tail only data above the 
median count rate are considered to compute variance and ex- 
pected variance. Note that the modified definition is still af- 
fected by the problem of "representativity" of all the source 
states. Since the median is always smaller than the average 
(for all energies, and for all trial time binnings), the values of 
F va r,med are larger than F var , but in every other respect the results 
are qualitatively unchanged. 

4.3.2. Doubling/halving timescale 

Similarly, the "minimum halving/doubling" timescale has 
been computed for each energy band taking the average of val- 
ues for the different input light curve binning (but discarding 
the 200 s one because too strongly subject to spurious results). 
As in Zhang et al. (1999) we rejected a T s horuy if the fractional 
uncertainty was larger than 20 %. 

In order to minimize the contamination by isolated data 
points we filtered the light curves excluding data lying at more 
than 3-<7 from the average computed over the 6 nearest neigh- 
bors (3 on each side), and for the 500 s binning light curves 
we did not consider pairs between data closer than 3 positions 
along the time series. Finally instead of taking the single abso- 
lute minimum value of T s horuy (over all possible i, j pairs) T s h or t 
is determined as the average over the 5 lowest values. 

The resulting T s h or t are listed at the bottom of Table 2 and 
shown in Figure 6a. The main findings are: 

1 . there is not significant difference between 1997 and 1998 
observations, at any energy. 

2. There is weak evidence that the softest X-ray band ex- 
hibits slower variability, while 

3. the timescales for the three higher energy bands are indis- 
tinguishable, and there is no sign of a trend with energy. 

We also distinguish doubling and halving timescales. Unfor- 
tunately it has not been possible to obtain an estimate of the 
doubling time for the lowest energy band for 1998 data, be- 
cause there are only a few, uncertain, data during the rise of the 
flare. We find that (see Table 2, and Figure 6b,c): 

4. again the softest X-ray band shows longer (halving) 
timescales; 

5. there is marginal evidence that the doubling timescale is 
shorter than the halving one (i.e. the rise is faster than the 
fall). 

This latter result seems to hold both for 1997 and 1998 light 
curves and for each single trial binning time, and for each en- 
ergy (except in 5 cases out of 35) 12 . However it is strongly 

12 We stress here that although the individual timescales are statistically con- 
sistent with each other, the behavior appear to be systematic at all energies. 
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Table 2 
Basic Variability Estimators 





LECS 0.1-0.5 


LECS 0.5-2 


MECS 2-3 


MECS 4-6 


Energy (keV) 


0.26 


1.18 


2.36 


4.76 




r.m.s. 


variability estimator F va r 




1998 All 


0.151 ±0.006 


0.267 ±0.009 


0.315 ±0.009 


0.366 ±0.011 


1997 All 


0.145 ±0.035 


0.215 ±0.007 


0.257 ±0.020 


0.305 ±0.012 


1997 sub-1 


0.154±0.039 


0.189±0.008 


0.198 ±0.014 


0.209 ± 0.024 


1997 sub-2 


0.176 ±0.010 


0.282 ±0.013 


0.340 ±0.012 


0.385 ±0.014 


1997 sub-3 


0.154 ±0.024 


0.289 ±0.018 


0.374 ±0.015 


0.432 ±0.017 


shortest timescale T s i, ort (ks) 


1997: all data 


80.6 ±46.6 


28.6 ± 6.4 


25.5 ±2.4 


25.5 ±2.0 


halving 


167.4±31.0 


62.9 ± 4.4 


31.0±3.1 


32.6 ±8.0 


doubling 


81.7±62.1 


28.2 ± 6.6 


27.3 ±4.1 


28.7 ±3.4 


1998: all data 


63.8 ± 12.1 


28.1 ± 4.2 


27.3 ±2.4 


21.4± 1.5 


halving 


63.8± 12.1 


35.7± 4.1 


36.9 ±8.7 


34.5 ±2.4 


doubling 


a 


30.4±12.1 


28.4 ±3.3 


23.6±1.9 



d It has not been possible to derive a reliable estimate for this case due to insufficient data. 



weakened by the averaging over the several time binnings be- 
cause T s hort slightly increases with the duration of the time bins, 
thus yielding a larger uncertainty on the average. 

5. 1998 DETAILED TEMPORAL ANALYSIS 

We performed a more detailed analysis of the characteristics 
of the April 21 st flare, in particular to determine any significant 
different behavior among different energy bands. Three are the 
main goals: i) to measure the flare exponential (or power-law) 
decay timescales (or power-law index); ii) to evidence possible 
differences in the flare rise and decay timescales; iii) to look for 
time lags among the different energies. 

We considered the energy ranges already adopted in Sec- 
tion 4. 

5.1. Decay Timescales 

A first estimate of the decay timescales has been presented 
by Maraschi et al. (1999), where a simple exponential law has 
been fit to the decaying phase of the light curves, and a depen- 
dence of the e-folding timescale on the energy band has been 
found. Here we analyze the issue of determining the timescale 
in more detail and in particular we consider the requirement for 
the presence of an underlying steady emission. 

In particular, the decay timescales for different energy bands 
have been estimated by fitting the post-flare light curves with 
an exponential decay superimposed to a constant (quasi-steady) 
flux eventually constrained by the fit (we immediately stress 
that an additional non-zero steady emission is necessary in or- 
der to obtain a meaningful fit with an exponential decay). In 
one set of fits the underlying contribution is set to zero. More 
precisely, we performed the fits with three possible conditions 
for the level of the steady contribution: un-constrained, con- 
strained and zero. 



The analytical expression of the function used in the fits is: 

; - T 1C f 

Fsteady + Fflaring <? T 



F(f) 



^steady 



l + TZe~ 



The model parameters are then the decay timescale r, the abso- 
lute value of F steac | y , and the ratio 1Z between the flaring and the 
steady component taken at a reference time T re f, which ideally 
should correspond to the peak of the flare. 

The "reference time" T re f for April 21 st flare was T re f = 25 x 
10 3 s (T= = 10924.0 TJD), and T ref = 170 x 10 3 s for April 
23 rd . Note that the value of 1Z obtained from the fit to the second 
dataset is only a lower limit on the actual flare amplitude since 
we do not have an estimate of the time at which the peak of this 
putative second flare occurred. 

In any case, as our focus is actually on the properties of the 
first, well defined, outburst, we used the parameters obtained 
from the April 23 rd data only to constrain the contribution from 
the steady component. Furthermore, it turns out that the deter- 
mination of the timescale r is not affected by the uncertainty on 
its contribution (see Table 3 and Fig. 8). 

We adopted a 1000 s binning for the 0.5-2 keV and 2-3 keV 
data, which have better statistics, while we the binning time for 
the 0. 1-0.5 keV and 4-6 keV light curves is of 2000 and 1500 s, 
respectively. The analysis has been performed in two stages, for 
each energy band: 

i) Fit to the April 21 st and 23 rd datasets independently, 
yielding values of r, F stea d y and 1Z. As an example, in 
Figure 7 the 2-3 keV light curve is shown together with 
its best fit models for April 2 1 st and 23 rd . 

ii) Fit of the April 21 st post-flare phase constraining the 
level of the underlying steady component to be below 
the level of April 23 rd (the actual constraint is the up- 
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FIG. 6. — (a) comparison of global T s h rt f° r 1997 (crossed squares) and 1998 
(diamonds), (b) comparison of doubling (black triangles) and halving (white 
triangles) timescales for 1997, and (c) comparison of doubling (black squares) 
and halving (white squares) timescales for 1998 (see text). 



per bracket for a two parameters 90% confidence inter- 
val). We assume that the underlying emission varies on a 
timescale longer than that spanned by the observations. 

As anticipated we also modeled the decay of the April 21 st flare 
setting to zero the offset. 

The results are summarized in Tables 3 and 4. The best fit 
r for all the three cases are plotted in Fig. 8 as a function of 
energy. 

The 0.1-0.5 keV timescale is not affected by the constraint 
on the baseline flux and the 4-6 keV confidence interval only 
suffers a minor cut (smaller errorbar), while both the 0.5-2 and 
2-3 keV parameters are significantly changed because the best 
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FIG. 7. — Example of the fit of the light curve with a constant+exponential 
decay. The 2-3 keV data are shown, with a 1000 s time bin, together with the 
best fit model (solid line through the data points), and the best fit level of the 
steady component (three dashed lines corresponding to the best fit value and 
the 90% confidence interval, for two interesting parameters). 



100 

,— , 90 
jo 80 

£ 70 



60 
50 



CO 

o 
w 

CD 

B 



40 - 



30 - 



20 



~| is 1 r*n — i i i i\j 1 1 1 — i — i i i i | 




Energy [keV] 

FIG. 8. — Decay timescales r plotted versus the barycentric energy of the 
data from which they have been derived. The different symbols are: crossed 
circles for zero-baseline, black triangles for the un-constrained baseline, and 
diamonds for the constrained baseline case. The dashed lines show the loci of 
the relation r ~ E~'/ 2 , that holds if pure radiative cooling is responsible for the 
spectral variability during the burst decay. 



fit un-constrained level of the steady contribution is signifi- 
cantly higher than allowed by the April 23 rd data (see for in- 
stance Fig. 7). 

The main findings are: 

1) the decay timescales depend critically on the presence/ 
absence of a contribution by a non variable component. 
This is true not only for the value themselves (for the 
case without baseline are between a factor 2 and 4 longer) 
but also for the relationship between the timescales and 
energy. In fact: 

la) In the cases with baseline, the timescales range be- 
tween 30 and 45 x 10 3 s, and do not show a clear (if 
any) relationship with the energy, rather suggesting 
an achromatic post-flare evolution. In fact, accord- 
ing to a x 2 tes t. the values of r for the 4 energies are 
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Table 3 

Parameters of fit to light curves^ 







LECS 0.1-0.5 


LECS 0.5-2 


MECS 2-3 


MECS 4-6 


constant + exponential decay 


April 21, zero-baseline 










r c (ks) 
April 21, un-constrained 


125.3+H 


91.8 + }J 


SI 0+1.8 


74.4+3;' 


Fsteady (Cts/s) 

r c (ks) 
April 21, constrained 


0.77*;* 
29.9+1; 4 . 


1.38+°;" 

9 90+O.I2 
z - -0.18 

34.4l;l 


1 co+0.22 
L5S -0.16 

1 01+ - 07 
32 1+ 43 


2 97+ 3 °l 
z - -0.77 

39+ 09 
u - -0.12 

40.7+1;' 


r steady 


(cts/s) 


77+<>- 24 
'-0.13 

0.67^ 07 
[< 0.72] 


l- 94 -o.07 
[< 1.94] 


7 IQ+0.16 

°- 7 %,03 
[< 0.79] 


2 97+ 3 °l 
z - -0.49 

0.39_ 012 
[< 0.39] 


r c (ks) 
April 23, un-constrained 


30 7+ 10 ' 6 


41 7 +L8 


43.0+20 


41.2+11 


Fsteady 

r c 


(cts/s) 
(ks) 


0.38+0;<» 
0.70*Jg 
29.0+^0 


95+ - 27 

u - -0.13 

1 77+O.I6 
1 -"-0.24 

66.6+211 


1 25+ - 24 
55.6+1 4 ; 


1 35+ ' 17 
-0.11 

36+ 03 

U - JO -0.03 

m ?+ 8 - 6 
Jy - z -6.0 


constant + power law decay 






> 1.715 


> 15.15 


> 23.04 


> 39.37 


r steady 


(cts/s) 


< 0.600 


< 0.619 


< 0.489 


< 0.175 






54+ 037 

"""-0.18 


54+ 005 

-0.02 


62+ 016 
-0.02 


70+ 09 
u -' u -0.03 



a In the un-constrained case the flux level of the steady component is left free to adjust during the fit. In the constrained case its 
flux level is required to be less or equal to the upper value of the 90 % confidence interval obtained by the fit of the April 23 rd light 
curve in the corresponding energy band. The values reported within square brackets represent the imposed upper limits on F ste ady 
(we do not quote the upper error value because in all four cases the upper boundary of the confidence interval simply coincides with 
the imposed upper limit). Quoted errors are for a l-o confidence interval for one or two interesting parameters (i.e. A\ 2 = 1.0 or 
2.3) for the case without and with the baseline, respectively. 

b 7 

We do not list the x values because they are not really meaningful estimates of the goodness of the fit. This happens because 
the errors on each bin of the light curves are very small and there is significant variability on very short timescales which of course is 
not accounted for by the simple exponential decay and eventually it has the effect of giving an artificially high \ 2 value. For the fits 
of April 21 st , the values of reduced \ 2 are between 1.6 and 4.1 for the cases including the baseline (for both exponential and power 
law decays), and between 2.6 and 6.5 for the pure exponential (zero-baseline) case. The number of degrees of freedom (for April 
21 st datasets) is between 20 and 35, depending on the energy band. 

c Timescale of the exponentially decaying flaring component. 

d Ratio between the flaring and steady components at T ref (see text). 

e 90 % confidence upper limit, for A\ 2 = 4.61, for the level of the steady component. 

power law decay spectral index. 
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Table 4 

tlmescales for different exponential decay fits 



LECS 0.1-0.5 LECS 0.5-2 MECS 2-3 MECS 4-6 



Tend = 


100 x 


10* 


Tend — 


90 x 


10 3 


Tend — 


80 x 


10 3 


Tend = 


70 x 


10 3 


Tend = 


60 x 


10 3 




pb 





April 21, zero-baseline 



125.3+£| 9L8!^ 83^™ 74.4+ 3 ,;' 



,+3.7 
-2.0 



120.3+ 8 ; 2 89.0+2;0 80.3+ 2 ; 4 , 72.0 

116.2+9;" 84.9+j;?, 73.8 +2 | 6S.5+ 3 ; 8 

106.2+7| 80.7+2° 70.8+|5 68 - 5 -3:o 

95.3+^ 76.0+|] 70.0+ 3 J 70.3+5;' 

4 x 10~ 2 6 x 10" 8 1 x 10" 5 



April 21, un-constrained baseline 



Tend — 


100 x 


10 3 


29.9+»; 4 


34 4+ 4 - 5 

-3.2 


32 1 +4 ' 3 


40.7+1 1 ;' 


Tend = 


90 x 


10 3 


33.8+1 4 ; 8 


37 7+ 5 - 5 


34.0+5;6 


46 3 + ' 6 - 5 


Tend — 


80 x 


10 3 


26.9+10;? 


34.5^;| 


38.6+1"; 4 




Tend = 


70 x 


10 3 


o-i -j+24.0 
-11.1 


34.6+9; ? 


45 6+ 26 - 5 




Tend = 


60 x 

pb 


10 3 


32 6 +m4 


26.3+J; , 


24.2+1'; 2 


18-9+1^ 



April 21, constrained baseline 



Tend = 


100 x 10 3 


30 7 +10 - 6 


41 7+ L8 


43.0! 2 ;? 


41.2 + 1°;« 


Tend — 


90 x 1() 3 


i-> o+15.1 
-8.0 


41 - 2 -L3 


42.0+ 3 


46.2 + ^;f 


Tend = 


80 x 10 3 


26.6+1'- 3 


40.8+" 


40.6+ 8 ;?, 


53 1+ 19 - 4 


Tend — 


70 x 1() 3 


3i-n.8 


40.5+6;' 


45.5+ 2 *; 2 


63 5 + 9 - 3 


Tend = 


60 x 10 3 

pb 


ii 7+69.3 
-9.4 


40.2+5;° 


41 5+ 9 - 6 


a a a+23.0 



a The decay of the April 21 st flare is modeled with an exponential function over a non- 
variable flux contribution. The listed numbers are the characteristic t's (expressed in ks) 
of the exponential decay, as a function of energy and T en d- Quoted errors are for a 1-cr 
confidence interval for one or two interesting parameters for the case without and with the 
baseline, respectively. The timescales for the T en d = 100 ks case are the same of Table 3. 

b Probability that the values of the timescale r obtained for the different choices of T cn d 
are drawn from a unique parent distribution (\ 2 with respect to the weighted average), for 
those cases where there is a significant rejection of the "null hypothesis". 
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consistent with coming from the same distribution 
(the so called "null hypothesis"). 

lb) On the contrary, in the case of pure exponential de- 
cay the timescales follow a weak inverse relation- 
ship with the energy, r oc E~° ls±0 02 (both when 
considering r or E as the dependent variable). The 
X 2 probability in favor of the "null hypothesis" is in 
this case negligible. 

2) the amplitude of the variability (as measured by 1Z) pos- 
itively correlates with the energy, as reported in many 
cases for HBLs (e.g. Sambruna et al. 1994, Ulrich et 
al. 1997). 

5.1.1. Is there a steady component ? 

Variability timescales and their energy dependence carry pre- 
cious information on the physics of the source: the puzzling de- 
pendence of the results on the hypothesis of the presence of a 
steady component requires further analysis. 

Here we only discuss a simple test of the robustness of the 
results of previous section, that however provides useful hints. 
This consists in determining r for different choices of the end 
time (T en d) of the interval over which the fit is performed. The 
underlying idea is that if the shape of the decay is close to the 
one we are assuming (exponential with/without an offset) one 
can expect that the model parameters do not change much as 
a function of T enc |. On the contrary, if the chosen model pro- 
vides only a poor description of the actual decay characteristics, 
the parameters will likely take different values depending on 
T en d, unless of course we postulate a process whose timescale 
changes with time. 

For each of the three cases (constrained, unconstrained, 
zero-baseline) we tried 5 different choices of the end time, 
namely T e „ d = 60, 70, 80, 90 and lOOx 10 3 seconds (from T ), 
the last one corresponding to the entire duration of the data cov- 
erage on April 21 st . 

To evaluate the likelihood that the 5 values of t obtained for 
these choices of T en d are drawn from the same (normal) dis- 
tribution, we considered the probability V of the \ 2 computed 
with respect to their weighted 13 average 14 . 

In Table 4 we summarize the results obtained under each of 
the different hypotheses. There are a few points to note: 

1. in the zero-baseline hypothesis the timescales obtained 
changing T en d are not consistent with being several dif- 
ferent realizations from a unique distribution. Further- 
more the clear trend of (decreasing) r with (decreasing) 
Tend supports the possibility that the variation of r is sig- 
nificant. 

2. On the contrary in both cases with baseline the values of 
t's for the 0.1-0.5, 0.5-2 and 2-3 keV bands are very 
close for all values of T en d- For the 4-6 keV band the 
scatter is somewhat larger, likely due to the fact that at 
this energy there is the largest amplitude of variability, 
i.e. the smallest contribution of the putative baseline. The 
same kind of effect of lowered sensitivity to the details of 
a putative baseline, is responsible for the good result of 
the 4-6 keV band in the case without baseline. 

13 To each measure we assigned as a (and in turn as weight = 1/cr 2 ), the 
mean between its + and - uncertainties. 

14 It should be noted that the weighted average has the property of minimizing 
the x 2 for a dataset, and thus the resulting values of probability that the data do 
belong to a unique distribution are upper limits. 




Energy [keV] 

FIG. 9. — Index of the power law decay plotted versus energy. Only the value 
for the April 21 post peak decay is plotted. 



3. Even in the case without baseline the values of r for 
the various energy bands become very similar for T en d = 
60 x 10 3 s. Therefore even if the failure of the test by 
the pure exponential decay model could be interpreted 
as the signature of a mechanism whose characteristic 
timescale changes with time (unlike a true exponential 
decay), there is anyway an indication that in the earlier 
stages of the decay the evolution is not energy dependent. 

On the basis of this simple consistency check we can therefore 
infer that if the decay were exponential there is evidence for the 
presence of an underlying non-variable component. Thus the 
post-flare phase does not show any dependence on energy, and 
the timescale is of the order of 30^40 ks. 

5.2. Power Law Decay 

As an alternative description of the decay of April 21 st flare 
we consider a power law, where Fflari n g oc (T/T le f)~ v , again leav- 
ing the possibility of an offset from a steady component. The 
results are summarized in Table 3 and Figure 9. 

A fundamental difference with respect to the exponential de- 
cay case is that an additional steady component is never re- 
quired: the power law fits yield only (not very stringent) up- 
per limits on F stea dy and only for the 0.1-0.5 keV band the for- 
mal best fit value for F s t ea dy is different from zero (F stea dy = 0.26 
cts/s). 

Again there is only a weak relationship of the decay prop- 
erties with energy, with a marginal indication that the flux de- 
crease is faster at higher energies. The values for the 4 energy 
bands are consistent with a common (weighted average) value 
of 77 = 0.583, while for the subset comprising only the three 
higher energy datasets the probability for the "null hypothesis" 
is reduced to about 0.07. 

5.3. Rise vs. Fall 

A further interesting question is whether the rise and the de- 
cay phases of the flare are characterized by the same timescale 
(i.e. if the flare is symmetric or asymmetric) and in particular 
how these properties might correlate with energy. 

The observational coverage of the rise is not good enough to 
apply the same direct technique used for the decay (§5.1), to 
constrain the timescale of an exponential rise, with all the pa- 
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Table 5 

Difference between Rise and Fall Timescales 1 * 



case b 


LECS 0.1-0.5 


LECS 0.5-2 


MECS 2-3 


MECS 4-6 


V c 


Fsteady — 

Fsteady t> est fit un-constrained 
Fsteady best fit constrained 


9 Q+33.4 
^■ y -12.3 

1 2+ 6 - 9 
-4.3 

1 3+ 6 ' 2 
1- -4.3 


26 3+ 15 
13 1 M 

"■'-1.9 
Mi.!** 


33.1 + »j 
14.8^-9 

19.8+W9 


38.6^j 
23.6*;* 
24.6!" 


3 x 10 -2 
1 x 10" 2 



a The values of AT = Tf a ii-T r i se have been computed for two different binning of the light curves, 500 and 1000 
s, and the average is shown (in ks). Confidence intervals are l-cr for 1 interesting parameter, because the baseline 
was not free to vary. The upper errorbars (giving the largest possible value of AT) derive from the upper boundary 
of the Tf a n confidence interval and the lower boundary of the T r ; sl , confidence interval. The lower errorbars (giving 
the lowest possible value of AT) derive from the lower boundary of the Tf a n confidence interval and the upper 
boundary of the T r j se confidence interval. 

b In each case the level of the steady component was kept fixed at the best fit value for each energy band, as 
listed in Table 3. 

c Probability that the four values of AT are drawn from the same parent distribution, computed using the \ 2 OI 
the data with respect to their weighted average. 



with fits to the light curves, are shown in Table 5, and plotted in 
Figure 10. The flare is symmetric for the softest energy band, 
while it is definitely asymmetric in the other three bands, with 
an asymmetry which appears to increase systematically with 
energy: at higher energies the rise phase is increasingly faster 
than the decay (the significance is reported in Table 5). 

It is worth stressing that these numbers indicate only relative 
changes in the timescales before and after the peak of the flare. 
On the other hand, since we found that the decay timescales are 
very similar once the baseline flux is taken into account (§5.1), 
this asymmetry probably means that the rise times get shorter 
with increasing energy. 



5.4. Time Lag 

We performed a detailed cross correlation analysis using two 
different techniques suited to unevenly sampled time series: the 
Discrete Correlation Function (DCF, Edelson & Krolik 1988) 
and the Modified Mean Deviation (MMD, Hufnagel & Breg- 
man 1992). Moreover, Monte Carlo simulations, taking into ac- 
count "flux randomization" (FR) and "random subset selection" 
(RSS) of the data series (see Peterson et al. 1998), were used to 
statistically determine the significance of the time lags between 
different X-ray bands obtained with the DCF and MMD. We 
refer for the relevant details of such analysis to Zhang et al. 
(1999). 

We binned the light curves over 300 s in the 0. 1-1 .5 and 3.5- 
10 keV bands (whose effective barycentric energies are E ~ 0.8 
and E~5 keV), and a trial time step of 720 s is adopted for 
both the DCF and MMD. DCF amplitude versus time lag is 
plotted in Figure 11. A negative value means that variations 
in the 3.5-10 keV band light curve lag those occurring in the 
0.1-1.5 keV one (i.e. hard lag). The best Gaussian fits for both 
DCF and MMD result in negative time lags of -2.8 ±0.2 (DCF) 
and -2.1 ±0.3 (MMD) ks, indicating that the medium energy 
X-ray photons lag the soft ones. 

The Cross Correlation Peak Distribution (CCPD) obtained 
from the FR/RSS Monte Carlo simulations is shown in Fig- 
ure 12 both for the DCF and the MMD methods. The average 
lags resulting from CCPD are -2.T_\l ks for DCF, and -2.3^ 
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FIG. 10. — Difference between the post-peak decay and the pre-peak rise 
timescales, plotted versus energy. The different symbols are: crossed circles 
for F stea< j y = case, black triangles for the F stell( j y set at the best fit value of the 
un-constrained case, and diamonds for F stea( ] y set at the best fit value obtained 
for the April 23 rd data (see Table 3). 



rameters left free. Thus the flux level of the steady component 
is fixed at the best fit value obtained for the decay, or set to 0. 

We are mainly interested in the comparison of the proper- 
ties of the outburst just before and after the top of the flare (see 
discussion in Paper 2), and the interval of about 17-18 ks pre- 
ceding the flare is long enough for this purpose. 

We thus focused on the "tip" of the flare, considering data 
on the decay side only up to the time when the flux reaches 
approximately the level at which the observation starts. The 
duration of the post-peak intervals needed to decrease the flux 
to the pre-flare level are of about 15-20, 25, 30-35 and 40 ks 
for the four energy bands. Since the uncertainties on the decay 
time of the light curves are not larger than 5 ks, the differences 
in duration and the trend (longer time at higher energy) are real. 

The differences between rise and fall timescales, estimated 
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FIG. 1 1 . — Result of the cross-correlation DCF method applied to the 0.1-1.5 
and 3.5-10 keV light curves. Negative values corresponds to the 3.5-10 keV 
band lagging the 0.1-1.5 keV one. The continuous curve superimposed to the 
data is the best fit with a constant+Gaussian model. The actual fit has been 
performed only to the "core" of the data, i.e. taking into account the data in the 
interval ranging from -30000 to 30000 seconds. 
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FIG. 12. — Cross Correlation Peak Distributions for the FR/RSS simulations. 
The solid line is for the analysis with DCF technique, and the dashed for the 
case of MMD. Negative values corresponds to the 3.5-10 keV band lagging the 
0.1-1.5 keV one. 



ks for MMD (1 a confidence intervals, two sided with respect 
to the average), confirming the significance of the above results 
with high (> 90 %) confidence. 

The total integral probabilities for a negative lag (that would 
be the actual measure of the confidence of the "discovery" of 
the hard lag) are ~ 95.0% (DCF) and ~ 98.7% (MMD). 

6. SUMMARY OF VARIABILITY PROPERTIES 

We presented a comprehensive temporal analysis of the flux 
variability characteristics in several energy bands, of BeppoSAX 
observations of Mkn 421. The primary results of this study are 
the following: 

1. The fractional r.m.s. variability is higher at higher ener- 
gies. 

2. The fractional r.m.s. variability does not change with the 
brightness state of the source. 

3. The minimum halving/doubling time is longer for the 



softest energy band. 

4. The minimum halving/doubling time at a given energy 
does not change with the brightness state of the source. 

5. There is a hint that the doubling timescale is shorter than 
the halving timescale. 

The findings 1 and 2, on F var , confirm the well known (although 
not fully explained) phenomenology of HBLs (see Ulrich et al. 
1997 and references therein). 

The results on the variability timescales (findings 3, 4, 5) 
show good agreement with those found with the more thorough 
analysis of the 1998 flare, that can be summarized as: 

6. The flare decay is consistent with being achromatic, both 
if modeled as an exponential decay with an additional 
contribution from a steady component and in the case of 
a power law decay. 

7. The results on timescales are at odds with the simplest 
possibility of interpreting the decay phase, i.e. that it is 
driven by the radiative cooling of the emitting electrons. 
This would give rise, in the simplest case, to a depen- 
dence of the timescale on energy, r ~ E" 1 / 2 . The tracks 
corresponding to this relationship are overlaid to the data 
in Fig. 8, and it is clear that they can not be reconciled. 

8. The harder X-ray photons lag the soft ones, with a delay 
of the order of few ks. This finding is opposite to what 
is normally found in the best monitored HBL where the 
soft X-rays lag the hard ones (e.g. Takahashi et al. 1996 
for Mkn 421; Urry et al. 1993, Zhang et al. 1999, and 
Kataoka et al. 2000 for PKS 2155-304; Kohmura et al. 
1994 for H0323+022). The latter behavior is usually in- 
terpreted in terms of cooling of the synchrotron emitting 
particles. 

9. Possible asymmetry in the rise/decay phases: the flare 
seems to be symmetric at the energies corresponding 
(roughly) to the peak of the synchrotron component (see 
the spectral analysis in Paper II), while it might have a 
faster rise at higher energies. This could be connected to 
the observed hard-lag. 

10. Finally, the characteristics of the April 21 st flare suggest 
the presence of a quasi-stationary emission contribution, 
which seems to be dominated by a highly variable peaked 
spectrum. 

7. CONCLUSIONS 

BeppoSAX has observed Mkn 421 in 1997 and 1998. We 
analyzed and interpreted the combined spectral and temporal 
evolution in the X-ray range. During these observations the 
source has shown a large variety of behaviors, providing us a 
great wealth of information, but at the same time revealing a 
richer than expected phenomenology. 

In this paper we have presented the first part of the analysis, 
focused on the study of the variability properties. 

The fact that the higher energy band lags the softer one (with 
a delay of the order of 2-3 ks) and the energy dependence of 
the shape of the light curve during the flare (with faster flare 
rise time at higher energies) provide strong constraints on any 
possible time dependent particle acceleration mechanism. In 
particular, if we are indeed observing the first direct signature 
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of the ongoing particle acceleration, progressively "pumping" 
electrons from lower to higher energies, the measure of the de- 
lay between the peaks of the light curves at the correspond- 
ing emitted frequencies would provide a tight constraint on the 
timescale of the acceleration process. 

The decomposition of the observed spectrum into two com- 
ponents (a quasi stationary one and a peaked, highly variable 
one) might allow us to determine the nature and modality of the 
energy dissipation in relativistic jets. 

In Paper II we complement these findings with those of the 
time resolved spectral analysis and develop a scenario to inter- 
pret the complex spectral and temporal phenomenology. 
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APPENDIX 

A. DEFINITIONS OF F VA r AND Tshort 

The fractional r.m.s. variability amplitude is a useful param- 
eter to characterize the variability in unevenly sampled light 
curves. It is defined as the square root of the so called excess 
variance (e.g., Nandra et al. 1997). This parameter, also known 
as the true variance (Done et al. 1992), is computed by taking 
the difference between the variance of the overall light curve 
and the variance due to measurement error, normalized by di- 
viding by the average squared flux (count rate). 

We consider a dataset Ffai) (i = l,N), with an uncertainty a,- 
assigned to each point. 

The fractional r.m.s. variability parameter is then defined as: 



F v ar — 



_ (4- A|) 1/2 

<*■> 



where 



a - 



1 N 



(Al) 

(A2) 
(A3) 
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The definition of the "shortest variability timescale" is the 
following: 



Tshort = min(T s h ort ,i/) 

ViJ 



1 short.;) — 



AF U 



(A4) 
(A5) 



where /vy = (fl+F,)/2, AF U = F-Fj, and AT U = 7} -7). 



